Neutrino-nucleus reactions and nuclear structure* 
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The methods used in the evaluation of the neutrino-nucleus cross section are reviewed. Results are 
shown for a variety of targets of practical importance. Many of the described reactions are accessible 
in future experiments with neutrino sources from the pion and muon decays at rest, which might 
be available at the neutron spallation facilities. Detailed comparison between the experimental and 
theoretical results would establish benchmarks needed for verification and/or parameter adjustment 
of the nuclear models. Having a reliable tool for such calculation is of great importance in a variety 
of applications, e.g. the neutrino oscillation studies, detection of supernova neutrinos, description 
of the neutrino transport in supernovae, and description of the r-process nucleosynthesis. 



I. INTRODUCTION 

There is now a convincing evidence that neutrinos are massive and the existence of neutrino oscillations have been 
convincingly demonstrated. This conclusion is based to a large extent on the observation of the neutrino induced 
reactions on complex nuclei, which also play essential roles in various aspects of nucleosynthesis as well as supernova 
collapse and supernova neutrino detection. Most of them have not been studied experimentally so far and their cross 
sections, which are needed in all applications, are typically based on nuclear theory. Spallation neutrino sources with 
their significant neutrino fluxes represent a unique opportunity to establish several benchmark measurements of the 
most significant neutrino-nucleus reactions. These measurements, in turn, can be used to gauge the accuracy and 
reliability of the corresponding nuclear models. 

In this work we review selected theoretical results of particular importance for neutrino detection, supernovae, and 
nucleosy thesis. Theoretical description of the neutrino induced reactions is a challenging proposition, since the energy 
scales of interest span a vast region, from the few MeV for solar neutrinos, to tens of MeV for the interpretation of 
experiments with the muon and pion decay at rest and the detection of supernova neutrinos, to hundreds of MeV 
or few GeV for the detection atmospheric neutrinos. While reactions induced by low-energy neutrinos are sensitive 
to details of nuclear structure, GeV neutrinos, like other weak probes of similar energy, interact dominantly with 
individual nucleons in the nucleus, which can then be treated as an ensemble of non- interacting but bound protons 
and neutrons. 

First, let us briefly review the general formalism adopted for the analysis of the charged current reactions, 



v e +zXjy -^z+i Xn_! 
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and its analogs with v e as well as with muon neutrinos and antineutrinos. The formalism can be easily modified for 
the neutral current reactions 
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In the derivation of the relevant cross sections we follow the prescription given by Walecka [1] which is based on 
the standard current-current form for the weak interaction Hamiltonian governing these reactions. After a multipole 
expansion of the weak nuclear current and application of the extreme relativistic limit (final lepton energy Eg » 
lepton mass mic 2 ) the neutrino (antineutrino) cross section for excitation of a discrete target state is given by [1,2]: 
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Here is the angle between the incoming and outgoing lcpton, and = (u), q) (q = \q\) is the four-momentum transfer. 
The minus-(plus) sign in Eq. (5) refers to the neutrino (antineutrino) cross section. The quantities Mj,Lj,jf and 
jmag j en0 |- e j- ne multipole operators for the charge, the longitudinal and the transverse electric and magnetic parts 
of the four-current, respectively. Following Refs. [1,3] they can be written in terms of one body operators in the 
nuclear many-body Hilbert space. The cross section involves the reduced matrix elements of these operators between 
the initial state Ji and the final state J/. (See Refs. [4,5] for the slightly more complicated formula valid also for 
nonrelativistic final lepton energy.) 

For low energy electrons and positrons the Fermi function F(Z, Eg) accounts for the Coulomb interaction between 
the final charged lepton and the residual nucleus in the charged-current processes. We use the Coulomb correction 
derived by numerical solution of the Dirac equation for an extended nuclear charge [6]: 



F(Z,E e ) = F {Z,E e ) ■ L , with F (Z,E e ) = A(2 Pl R) 2 ^ 
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Here Z denotes the atomic number of the residual nucleus in the final channel, Eg the total lepton energy (in units of 
mgc 2 ) and pi the lepton momentum (in units of mgc), R is the nuclear radius (in units of -^—) and 7 and y are given 

* ■ ' ■ Til p_ C 

by (a = fine structure constant): 
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The numerical factor Lq in (6), which describes the finite charge distribution and screening corrections, is nearly 
constant(« 1.0), and can be well approximated by a weakly decreasing linear function in pg. 

At higher energies, and for muons at essentially all energies, the Fermi function valid for s-wave leptons is a poor 
approximation for the Coulomb effect since higher partial waves also contribute for pgR > 1. Guided by the distorted- 
wave approximation of quasielastic electron scattering, we treat in that case the Coulomb effects in the 'Effective 
Momentum Approximation' in which the outgoing lepton momentum pg is replaced by the effective momentum 
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where Vc(0) = 3e 2 Z/2R is the Coulomb potential at the origin. In the work presented here, the Coulomb effect is 
taken into account not only by using the effective momentum, but also by replacing the phase space factor pgEg by 
p e ff E e ff (see also [7] where this procedure is called modified effective momentum approximation, and shown to work 
quite well). In practice we use a smooth interpolation between these two regimes of treatment of the Coulomb effects. 

We calculate the differential cross section (3) as a function of the initial neutrino energy e„, the excitation energy 
of the nucleus u> and the scattering angle 6. The three-momentum transfer q = \q\ is equal to 



(E v - pg) 2 + 4E„ ■ pg sin 2 ® 



lu 2 + 4£„ • {E v - u) sin 2 ® , 



(9) 



where the last expression is valid in the relativistic limit (Eg 3> mg) for the final lepton. The total cross section is 
obtained from the differential cross sections by summing (or integrating) over all possible final nuclear states and by 
numerical integration over the angles. 

As explicitly used in the derivation of the cross section formula above, in neutrino-induced reactions the nucleus is 
excited by multipole operators 0\ which scale like (qR/hc) x , where R is the nuclear radius (R ~ 1.2A 1 / 3 fm). As the 
momentum transfer is of the order of the neutrino energy E v , neutrino- nucleus reactions involve multipole operators 
with successively higher rank A with increasing neutrino energy. Since the nuclear Hamiltonian does not commute 
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with Ox, the response of the operator is fragmented over many nuclear states. However, for each multipole most of 
the strength resides in a collective excitation, the giant resonance, with a width of a few MeV. The centroids of the 
giant resonances of the various multipoles grow in energy with increasing rank, roughly like Xhtv, where Tiuj ~ 41 /A 1 / 3 
MeV is a typical energy shell splitting in the nucleus. Furthermore, since the phase space is proportional to piEg the 
higher outgoing lepton energies are preferred. This suggests that the average nuclear excitation energy to lags behind 
the increasing neutrino energy, i.e., for sufficcntly large neutrino energies the inequality u> < E v ~ Et holds. As a 
consequence, for neutrino energies significantly larger than the energies of the corresponding giant resonances, the 
neutrino-induced cross sections will depend on the total strength of the multipole excitation and its centroid energy, 
but will be less sensitive to its detailed energy distribution. Finally, at very high neutrino energies the neutrino will 
see the nucleus as an ensemble of independent bound nucleons and will interact with individual nuclcons. In that 
regime we will consider only the quasielastic channel in which the struck nucleon is ejected. This is the channel most 
widely used in neutrino detectors. 

From these general considerations we can identify three different energy ranges with quite different demands on the 
details with which the nuclear structure should be treated: 

i) For relatively low neutrino energies, comparable with the nuclear excitation energy, neutrino-nucleus reactions 
are very sensitive to the appropriate description of the nuclear response. Thus, low energy neutrino scattering requires 
a nuclear model which reproduces the important correlations among nucleons. The model of choice is the nuclear 
shell model, which accounts for nucleon-nucleon correlations via an effective interaction within a fixed model space 
for the valence nucleons. Nowadays, complete diagonalization for the lowest states in medium-mass nuclei (up to 
A <~ 60) is achievable in complete OTiui model spaces, i.e., considering all configurations of the valence nucleons in a 
full harmonic oscillator shell. For lighter nuclei, like 16 O, complete diagonalization can be performed in larger model 
spaces (ATwj for 16 0, 6Tioj for 12 C). Importantly, the shell model calculations have been proven to indeed reproduce the 
allowed (Fermi, Gamow- Teller) response for those nuclei for which diagonalizations can be performed in sufficiently 
large model spaces. For the lighter nuclei, where multi-shell calculations can be performed, the shell model also nicely 
describes forbidden transitions. An overview of recent shell model developments and applications is given in [8]. 

ii) The Random Phase Approximation (RPA) has been developed to describe the collective excitation of a nucleus 
by considering the one-particle one-hole excitations of the correlated ground state. In the standard RPA, all excited 
states are treated as bound states, leading to a discrete excitation spectrum. In the Continuum RPA (CRPA)(see 
e.g. [9,10] and references therein) the final states have the appropriate scattering asymptotics for energies above the 
nucleon-emission thresholds; consequently the excitation spectrum in the CRPA is continuous. The RPA or CRPA 
are the methods of choice at intermediate neutrino energies where the neutrino reactions are sensitive dominantly to 
the total strength and the energy centroids of the giant resonances. 

iii) At high incoming energies neutrinos, like other electroweak probes, scatter 'quasi-freely' on individual nucleons. 
The remaining nucleons can be treated as (non-interacting) spectators. This situation is realized in the Fermi gas 
model [11] where a full relativistic treatment of the hadronic weak vertex is included. The Fermi motion and the 
binding energy are characterized by parameters that can be determined from electron scattering experiments. In 
this approximation the expansion in multipoles is no longer necessary. The nuclear form factors for the quasi-free 
scattering are expressed in terms of the single nucleon matrix elements which depend only on the four-momentum 
transfer q 2 and on the nuclear momentum distribution. 

Above we pointed out the important role of collective excitations. The centroid position of these excitations deviates 
noticeably from the independent particle estimate of Xhu> due to the residual particle-hole interaction. Thus, empirical 
evidence or nuclear structure calculation is needed to determine their energy On the other hand, the total strength 
is often fixed by sum rules in an essentially model independent way. Well known examples of such sum rules is the 
Ikcda sum rule for the Gamow- Teller strength 

B(GT; Z -> Z + 1), - B(GT; Z -► Z - 1) 4 = 3{N - Z) , (10) 

i i 

or the Thomas-Reiche-Kuhn sum rule for the dipole strength 

i ^ 

Even though both of these sum rules could be violated to some extent (e.g., when the internal structure of the nucleons 
is not properly treated), the dependence on the neutron and proton numbers N and Z remains valid. Note that the 
Ikeda sum rule involves the difference of the strengths. However, in nuclei with neutron excess, N > Z, the second 
term, the total strength B(GT; Z — > Z — 1) is much smaller than the first one, and so the sum rule determines the 
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strength in the (p, n) channel. In symmetric nuclei with N = Z the Ikeda sum rule, unfortunately, does not help in 
fixing the Gamow- Teller strength. 

To calculate the various partial neutrino-induced reaction cross sections for neutrino-induced reactions we assume 
a two-step process. In the first step we calculate the charged current (V;, l~) and {pi, l + ) cross sections (where I — e or 
/i), or the neutral current cross section (u, v') as a function of excitation energy in the final nucleus. These calculations 
are performed within the RPA or CRPA and considering all multipole operators up to a certain J and both parities. 
In the second step one calculates for each final state with well-defined energy the branching ratios into the various 
decay channels using the statistical model code SMOKER [12]. As possible final states in the residual nucleus the 
SMOKER code considers the experimentally known levels supplemented at higher energies by an appropriate level 
density formula [12]. Proton, neutron, a and 7 emission are included in the code as decay channels. If the decay 
leads to an excited level of the residual nucleus, the branching ratios for the decay of this state is calculated in an 
analogous fashion [13]. Keeping track of the energies of the ejected particles and photons during the cascade, and 
weighting them with appropriate branching ratios and the corresponding primary charged- or neutral-current cross 
sections, we determine the various partial particle emission cross sections. 

II. COMPARISON OF DIFFERENT METHODS 

Here we demonstrate, using the neutrino interaction with 16 as an illustration, how different theoretical methods 
can be used at different neutrino energies. We show that at certain transition energy intervals the corresponding 
methods give essentially identical results. 

A. Shell model versus CRPA 
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FIG. 1. Comparison of the CRPA (full lines) and shell model (dashed lines) cross sections. The upper panel is for the f e 
induced reaction and the lower one is for the reaction induced by v e . 
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We first consider the description of the charged current reactions on 16 at relatively low energies. The reaction 
thresholds are 15.4 MeV for the 16 0(> e , e~) 16 F (which is unbound) and 11.4 MeV for the 16 0(9 e , e+) 16 N reaction. 

The shell model evaluation of the cross sections was performed back in 1987 by Haxton [14]. In that work the low 
lying positive parity states were described in a full 2 Two shell model. The transitions to negative parity states were 
described using the effective density matrices, scaled to describe measured form factors from electron scattering. 

The shell model results can be compared to the CRPA. The CRPA calculations used the finite range residual force 
based on the Bonn potential, and all multipole operators with J < 9 and both parities were included. The free 
nucleon form factors were used, with no quenching. The procedure was tested by evaluating the total muon capture 
rates (dominated by the negative parity multipoles) for 12 C, 16 O and 40 Ca [20], as well as the partial capture rates to 
the bound CP, 1~ and the 2~ ground state in 16 N [20,16]. Good agreement with these muon capture rates tests the 
method at momentum transfer q ~ m M ^100 MeV. 

We compare the cross sections evaluated by the two methods in Fig. 1 , where we show the cross sections evaluated 
in both methods and averaged over the Fermi-Dirac distribution corresponding to the temperature T and vanishing 
chemical potential. The agreement is excellent suggesting that both methods are capable of describing the weak 
reaction rates in this energy regime, provided that they can be successfully tested on relevant quantities, such as 
the muon capture rates, nuclear photoabsorption cross section, or inelastic electron scattering leading to the states 
populated by the weak processes. 

The angular distribution of the emitted electrons with respect to the incoming neutrino beam is shown in Fig. 2. 
Note the electron emission is predominantly in the backward direction at low energies (also obtained in the nuclear 
shell model), but it gradually changes to the forward one at higher energies. Thus, for E v > 500 MeV the direction 
of the electron can be used to determine the direction of the incoming neutrino. 
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FIG. 2. The CRPA angular distributions. 



We have thus identified the energy region, somewhere near about 50 MeV of neutrino energy, where the two discussed 
methods, the nuclear shell model and the CRPA, give essentially identical results. For lower energies the nuclear shell 
model is the method of choice. As the energy increases, the shell model calculations become increasingly difficult. 
The number of states increases rapidly, and the effective interaction to be used becomes more uncertain. However, 
as we argued above, at higher neutrino energies, above, say, E v > 100 MeV, the details of the nuclear correlations 
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become less important and what matters are the positions and strengths of giant resonances. The CRPA is capable 
of describing these quantities and thus, in our opinion, it is the method of choice at the intermediate neutrino energy 
range, approximately 100 MeV < E v < 500 MeV. As the energy increases further, the CRPA calculations become 
computationally more difficult (more multipoles and higher nuclear excitation energies must be included). At the 
same time, the nuclear response, at least for the quasielastic regime, becomes simpler. 

B. CRPA versus Relativistic Fermi Gas model 

The CRPA and Relativistic Fermi Gas model (RFG) methods are compared in Fig. 3 and agree remarkably well in 
both the total quasielastic cross section, as well as in the angular distribution of the outgoing electrons. The latter is 
particularly important, because the zenith angle distribution of the atmospheric neutrinos is based on the assumption 
that one can deduce the incoming neutrino direction (and hence its flightpath) from the direction of the observed 
charged lepton. Our CRPA calculations confirm that, indeed, below about 500 MeV of neutrino energy the emitted 
electron (or muon) are essentially uncorrelated with the direction of the incoming neutrino. Above these energies, the 
emitted lepton moves dominantly in the direction of the incoming neutrino, hence one can, statistically, correlate the 
two. This tendency, naturally, becomes more pronounced at higher E v values. 
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FIG. 3. Comparison of the relativistic Fermi gas model (full lines) and the CRPA (dashed lines). The parameters of the 
Fermi gas model were Pf — 225 MeV and e;, = 27 MeV. The cross sections are shown in the upper panel for the two indicated 
reactions. The angular distributions for E v = 300 and 500 MeV are shown in the lower panel. The CRPA results are shown 
with full lines, RFG with dashed lines. 

Let us stress that no attempt was made to adjust the Fermi gas model parameters to obtain the agreement 
demonstrated in Fig. 3; these are just standard values of the Fermi momentum pf and of the parameter e?, which 
characterizes the average nucleon binding energy. 
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Thus, we conclude that at energies E v ~ 300 MeV one can safely switch from the CRPA to the simpler Fermi gas 
model description of the quasielastic charged current reactions. However, it should be stresses that at low energies 
there are important checks in the form of muon capture rates as well as (for 12 C) neutrino induced reaction for the 
neutrino beams from the muon and pion decay at rest. Similar tests do not exist, or were at least less exploited, for 
~ 1 GeV neutrino energies. 

III. NEUTRINO 12 C INTERACTION 

The nucleus 12 C is particularly important for the study of the neutrino- nucleus scattering. Liquid scintillator 
detectors, for example KARMEN and LSND, contain hydrogen and 12 C nuclei. Therefore, a number of experimental 
results exist in this case as byproducts of the neutrino oscillation searches performed with these detectors. 

The measurements include charged-current reactions induced by both electron- [17,18] and muon-neutrinos [18], 
exciting both the ground and continuum states in 12 N. As discussed below, the inclusive cross section for 12 C(V e , e) 12 N* 
with the v e from the muon decay at rest (DAR) [17-19], agrees well with calculations, while in contrast, there is a 
discrepancy between calculations [20,21,23,24] (with some notable exceptions [25,26]) and the measured [18] inclusive 
cross section for 12 C(f M , /i) 12 N*, which uses higher energy neutrinos from pion decay-in-flight (DIF). The disagreement 
is somewhat disturbing in light of the simplicity of the reaction and in view of the fact that parameter- free calculations, 
such as those in [20-22] , describe well other weak processes governed by the same weak current nuclear matrix elements. 
Moreover, as shown in the following subsection, the exclusive reactions populating the ground state of the final nucleus, 
12 C(i/ e , e) 12 N ffs and 12 C(i^, ^) 12 N gs , and the neutral current reaction 12 C(^ e , v' e ) 12 C(15.11 MeV) have been measured 
[17,18] as well, and agree well with theoretical expectations. 

A. Exclusive reactions 

Among the states in the final nucleus 12 N, which is populated by the charged current reactions with beams of v e 
or z^j, the ground state V = 1 + plays a special role. It is the only bound state in 12 N, and can be recognized by its 
positron decay (T1/2 = 11 ms) back to 12 C. Moreover, the analog of the 12 N 9;S , the I* = 1+ state with isospin T = 1 
at 15.11 MeV in 12 C, can be populated by the neutral current neutrino scattering, and is recognizable by its emission 
of the 15.11 MeV photon. Finally, even though there are several bound states in 12 B, its ground state, the analog 
of the other two (I n ,T) = (1 + ,1) states, is the state most strongly populated in muon capture on 12 C. Again, the 
population of the bound states in 12 B can be separated from the continuum by observing its electron decay (T\/ 2 = 
20.2 ms). 

Theoretical evaluation of the exclusive cross sections is constrained by the obvious requirement that the same 
method, and the same parameters, must also describe the related processes, i.e. the positron decay of 12 N, the j3 
decay of 12 B, the Ml strength of the 15.11 MeV state in 12 C, and the partial muon capture rate leading to the 
ground state of 12 B. It turns out that this requirement essentially determines the neutrino induced cross section for 
the energies of present interest. It does not matter which method of calculation is used, as long as the constraints are 
obeyed. 

The comparison between the measured and calculated values is shown in Table I. There, three rather different 
methods of calculation were used, all giving excellent agreement with the data. 

The first approach is a restricted shell-model calculation. Assuming that all structure in the considered low-lying 
states is generated by the valence nucleons in the p-shell, and that the two-body currents (pion-exchange currents) are 
negligible, there are only four one-body densities (OBD) which fully describe all necessary nuclear matrix elements. 
In this case, it is necessary to use the one-body densities chosen (ad hoc) in such a way that all the auxiliary data 
mentioned above are correctly reproduced. This then gives the results listed in line 4 of Table I. 

Effects of configurations beyond the p shell might explain the need for the renormalization of the one-body densities 
produced by a reasonable p-shell Hamiltonian. Therefore, the rates of all the reactions are also evaluated in the 
Random Phase Approximation (RPA), which does include multishell correlations, while treating the configuration 
mixing within the p shell only crudely. Again an adjustment is needed (a "quenching" of all matrix elements by 
an universal, but substantial, factor 0.515). However, the neutrino cross sections in line 5 of Table I agree with the 
measurements very well. 

The third approach is the "elementary-particle treatment" (EPT). Instead of describing nuclei in terms of nucleons, 
the EPT considers them elementary and describes transition matrix elements in terms of nuclear form factors deduced 
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TABLE I. Comparison of calculated and measured cross sections, in units of 10~ 42 cm~ 2 and averaged over the corresponding 
neutrino spectra, for the neutrino induced transitions 12 C 9S — > 12 Ng 3 and 12 C 9S ^ 12 C(15.11 MeV). For the decay at rest the f e 
spectrum is normalized from E v = 0, while for the decay in flight the and spectra are normalized from the corresponding 
threshold. See the text for explanations. 





12 C(^ e ,e-) 12 N 9S 


12 C( W -) 12 N 9S 


12 C(z/, */) 12 C(15.11) 




decay at rest 


decay in flight 


decay at rest 


experiment [17] 


9.4±0.5±0.8 




11±0.85±1.0 


experiment [18] 


9.1±0.4 ± 0.9 


66±10± 10 




experiment [19] 


10.5±1.0 ± 1.0 






Shell model [27] 


9.1 


63.5 


9.8 


CRPA [20,21] 


8.9 


63.0 


10.5 


EPT [28] 


9.2 


59 


9.9 




1 J 1 ' 1 ' ' ' 1 

100 140 180 220 

Neutrino energy (MeV) 

FIG. 4. Energy dependence of the cross section for the reaction 12 C + -^ 12 N 9 . S + pT . 



from experimental data. The EPT approach was extended in Ref. [27] to the higher neutrino energies relevant to the 
LSND decay- in-flight f M 's by appropriately including the lepton mass. 

An example of the energy dependence of the exclusive cross section is shown in Fig. 4 for the f M induced exclusive 
reaction. As one can see, the cross section raises sharply from its threshold (E t hr = 123 MeV) and soon reaches its 
saturation value, i.e., it becomes almost energy independent. This means that the yield of the 12 C + reaction 
essentially measures just the flux normalization above the reaction threshold. At the same time, the yield is insensitive 
to the energy distribution of the muon neutrinos in the beam. 



B. Inclusive reactions 



The inclusive reactions 12 C(V e , e) 12 N*, with v e neutrinos from the muon decay-at-rest and Ciy^, /i) 12 N* with 
the higher energy neutrinos from the pion decay-in-flight populate not only the ground state of 12 N but also the 
continuum states. The corresponding cross sections involve folding over the incoming neutrino spectra and integrating 
over the excitation energies in the final nucleus. By convention, we shall use the term "inclusive" for the cross section 
populating only the continuum (i.e., without the exclusive channel) for 12 G{v e , e) 12 N* with the decay-at-rest v e , while 
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370 860 250 240 230 220 210 200 190 180 170 
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u [MeV] 

FIG. 5. Data (points with error bars) and calculated cross section for the inclusive electron scattering on 12 C as a function 
of the excitation energy u>. The corresponding momentum transfer is displayed on the upper scale. 



for the reaction 12 C(^ M ,^) 12 N* with the decay-in-flight the term is used for the total cross section (the exclusive 
channel then represents only a small fraction of the total). 

Muon capture, 12 C(/z, i/ M ) 12 B*, belongs also to this category. It involves momentum transfer of q « m^, intermediate 
between the two neutrino capture reactions above. Since 12 B and 12 N are mirror nuclei, all three reactions should be 
considered together. In this case again the term "inclusive" will be used only for the part of the rate populating the 
continuum in 12 B. 

Which theoretical approach should one use in order to describe such reactions? One possibility is to use the 
continuum random phase approximation (CRPA). The method has been used successfully in the evaluation of the 
nuclear response to weak and electromagnetic probes [9]. In particular, it was tested, with good agreement, in the 
calculation of the inelastic electron scattering [29] on 12 C involving very similar excitation energies and momentum 
transfers as the weak processes of interest. As an example Fig. 5 shows the comparison of the experimental data 
and the results of the CRPA for the inclusive electron scattering [30]. One can see that the CRPA describes quite 
well both the magnitude and shape of this cross section over the entire range of excitation energies and momentum 
transfers. 

For muon capture the CRPA [20] gives the inclusive rates of 0.342, 0.969, and 26.2 xlO 5 s" 1 for 12 C, 16 and 40 Ca; 
to be compared with the measured rates of 0.320, 0.924, and 25.6 xlO 5 s^ 1 for the same nuclei. This good agreement 
is again obtained without any parameter adjustment. In particular, as discussed in Ref. [20], no renormalization of 
the axial vector coupling constant <?a in nuclear medium is required. 

What are the momentum transfers and excitation energies involved in the inclusive reactions which we would like 
to describe? For the 12 C(^ e , e) 12 N* with the electron neutrinos originating in the muon decay at rest, the typical 
momentum transfer is (\q\) ~ 50 MeV, and the typical excitation energy is uj ~ 20 MeV. For the inclusive muon 
capture 12 C(/z~, ^) 12 B* we have (|<f|) ~ 90 MeV and the typical excitation energy is w ~ 25 MeV. Finally for the 
12 C(z^, /i - ) 12 N* with the muon neutrinos originating in the pion decay in flight at LAMPF we have (|<f|) ~ 200 MeV 
and the typical excitation energy is uj ~ 40 MeV. The excitation energies should be compared with the nuclear shell 
spacing fiuo ~ 41/A 1 / 3 MeV, which for 12 C is equal to about 18 MeV. Thus, in order to describe all the above inclusive 
processes in the framework of the nuclear shell model, one would have to include fully and consistently at least all 
2%ui excitations, and possibly even the 3Hoj ones. This is not impossible, but represents a formidable task. On the 
other hand, the CRPA can easily handle such configuration spaces. Moreover, it properly describes the continuum 
nature of the final nucleus. Finally, as argued above, the crudeness with which the correlations of the p shell nucleons 
is treated in the CRPA is expected to be relatively unimportant. 

For the inclusive reaction 12 C(z^ e , e~) 12 N*, with v e neutrinos from the muon decay-at-rest the calculation gives 
[20] the cross section of 6.3 x 10~ 42 cm 2 using the Bonn potential based G-matrix as the residual interaction, and 
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5.9 xl0~ 42 cm 2 with the schematic Migdal force. (The two different residual interactions are used so that one can 
estimate the uncertainty associated with this aspect of the problem.) Both are clearly compatible with the measured 
values of 6.4 ± 1.45[stat] ± \A[syst] x 10~ 42 cm 2 by the Karmen collaboration [17] (the more recent result gives 
somewhat smaller value 5.1 ± 0.6 ± 0.5 [50]) and with 5.7 ± 0.6[stat] ± 0.6[syst] x 10~ 42 cm 2 obtained by the LSND 
collaboration [18] . If one wants to disregard the error bars (naturally, one should not do that), one can average the 
two calculated values as well as the two most recent measurements and perhaps conclude that the CRPA calculation 
seems to exceed the measured values by about 10-15%. A similar tendency can be found, again with some degree of 
imagination, in the comparison of the muon capture rates discussed earlier. 

So far we have found that the CRPA describes the inclusive reactions quite well. Other theoretical calculations, 
e.g. [24,26] describe these reactions with equal success. This is no longer the case when we consider the reaction 
12 C(^ M , /i) 12 N* with the higher energy neutrinos from the pion decay-in-flight. This reaction involves larger 
momentum transfers and populates states higher up in the continuum of 12 N. The CRPA calculation [20,21] gives the 
cross section of 19.2 xlO -40 cm 2 , considerably larger than the measured [18] value of 11.3 ± 0.3[siai] ± l.8[syst] in 
the same units. The origin of the discrepancy is not clear, but as stressed in the discussion of the exclusive reaction, 
the flux normalization is not a likely culprit. While Ref. [23] confirms our result, Ref. [25] gets a value close to the 
experiment by using a generalization of the EPT approach. 

Other recent theoretical calculations span the region between the CRPA and experiment. So, Singh et al. [24] give 
16.65 ± 1.37 xlO -40 cm 2 , clearly higher than the experiment but somewhat lower than the CRPA. On the other hand, 
Ref. [26] gives 13.5 - 15.2 in the same units, a value which is even closer to the experiment. The main difference in 
that work is the inclusion of pairing which is not expected to represent a substantial effect. 

This discrepancy has been with us for quite some time now. It clearly exceeds the 10-15% perhaps suggested by 
the lower energy inclusive reactions discussed above. It would be very important to perform a large scale shell model 
calculation, including up to 3hu> excitations, to put the matter to rest. Attempts to do that are in Refs. [31,32]. 



IV. SUPERNOVA NEUTRINOS 



One of the most important application of the neutrino- nucleus interaction is the detection of supernova neutrinos. 
In this section, after few introductory remarks, we describe several examples of the calculated charged and neutral 
current cross sections on oxygen, argon, iron, and lead. All these nuclei are being considered (or actually are already 
used) as targets for the supernova neutrino detection. Some of these cross sections are amenable to tests using the 
spallation neutron source, since the neutrino spectra of stopped pions and muons are quite similar to the expected 
neutrino spectra from the core collapse supernovae. The general review of the field can be found in Ref. [33]. 

Supernova neutrinos from SN1987a, presumably all D e , had been observed by the Kamiokande and 1MB detectors 
[34,35] and have confirmed the general supernova picture. However, the supernova models predict distinct differences 
in the neutrino distributions for the various flavors and thus a more restrictive test of the current supernova theory 
requires the abilities of neutrino spectroscopy by the neutrino detectors. Current (e.g. Superkamiokande, SNO, 
KamLAND) and future detectors (including the proposed OMNIS [36] and LAND [37] projects) have this capability 
and will be able to distinguish between the different neutrino flavors and determine their individual spectra. For 
the water Cerenkov detectors (SNO and Superkamiokande) v x neutrinos can be detected by specific neutral-current 
events [13,38], while the OMNIS and LAND detectors are proposed to detect neutrons spalled from target nuclei by 
charged- and neutral-current neutrino interactions. 

Theoretical models predict characteristic differences in the neutrino distributions for the various neutrino flavors 
(so-called temperature hierarchy). The /i and r neutrinos and their antiparticles (combined referred to as v x ) decouple 
deepest in the star, i.e. at the highest temperature, and have an average energy of E v ~ 25 MeV. The v e and v e 
neutrinos interact with the neutron- rich matter via v e + n — > p + e~ and v e + p — > n + e + ; the D e neutrinos have a 
higher average energy (E v ~ 16 MeV) than the v e neutrinos (E v ~ 11 MeV). Clearly an observational verification 
of this temperature hierarchy would establish a strong test of the supernova models. The distribution of the various 
supernova neutrino species is usually described by the pinched Fermi-Dirac spectrum 

1 E 2 

n{Ev) = F 2 {a)T3 exp{(E„/T)-a} + l (12) 

where T, a are parameters fitted to numerical spectra, and F2 (a) normalizes the spectrum to unit flux. The transport 
calculations of Janka [40] yield spectra with a ~ 3 for all neutrino species. While this choice also gives good fits to 
the v e and D e spectra calculated by Wilson and Mayle [41], their v x spectra favor a = 0. In the following we will 



10 



present results for charged- and neutral current reactions on several target nuclei for both values of a. In particular 
we will include results for those (T, a) values which are currently favored for the various neutrino types (T in MeV): 
(T,a) — (4,0) and (3,3) for v e neutrinos, (5,0) and (4,3) for v e neutrinos and (8,0) and (6.26,3) for v x neutrinos. 
However, it is worthwhile to point out that the degree of separation in energy of the different flavors is somewhat 
model dependent, as shown e.g. in Ref. [39]. It is therefore even more important to determine the relevant parameters 
experimentally. 

As stated above, it is usually sufficient to evaluate the various neutrino-induced reaction cross sections within 
the RPA. However, the RPA often does not recover sufficient nucleon-nucleon correlations to reliably reproduce the 
quenching and fragmentation of the Gamow- Teller (GT) strength distribution in nuclei. Therefore, the response of 
the W = 1 + operator should be evaluated on the basis of an interacting shell model, if such calculations are feasible. 
While the double-magic nucleus 16 does not allow GT excitations, shell model calculations in reliably large model 
spaces are possible for 40 Ar and 56 Fe. For 208 Pb GT transitions are Pauli-blocked for (D e , e + ) reactions. The modeling 
of GT transitions in the (v e , e~) reactions on 208 Pb would require much too large model spaces; these transitions must 
also be evaluated within the RPA approach. 

In the following we will refer to a 'hybrid model' if the allowed transitions have been studied based on the interacting 
shell model, while the forbidden transitions were calculated within the random phase approximation. The studies for 
40 Ar and 56 Fe are performed in such a hybrid model, while the ones for 16 and 208 Pb are RPA calculations for all 
multipoles. 

We note that the GT operator corresponds to the A T = 1 + operator only in the limit of momentum transfer q — > 0. 
As it has been pointed out in [42,43], the consideration of the finite-momentum transfer in the operator results in 
a reduction of the cross sections. To account for the effect of the finite momentum transfer we performed RPA 
calculations for the X n = 1 + multipole operator at finite momentum transfer q (i.e. X(q)) and for q = (i.e., X(q = 0)) 
and scaled the shell model GT strength distribution by the ratio of X(q) and X(q = 0) RPA cross sections. The 
correction is rather small for v e neutrinos stemming from muon-decay-at-rest neutrinos (e.g., for LSND and Karmen) 
or for supernova v e neutrinos. The correction is, however, sizeable for higher neutrino energies. 

A. Oxygen 

Observation of neutrinos from the SN1987A did not allow to test in detail the neutrino distribution and, in particular, 
it gave no information about v x neutrinos which, as we discussed above, should decouple deepest in the star. The 
observability of supernova neutrinos has significantly improved since the Supcrkamiokande (SK) detector, with a 
threshold of 5 MeV and with 30 times the size of Kamiokande, became operational [44]. 

Clearly, many of the primary neutral- and charged-current i/-induced reactions in SK occur on 16 O. Atmospheric 
neutrinos and also supernova v x neutrinos have high enough energies so that the final nucleus in the primary reaction 
will be in an excited state which will then decay by a cascade of particle and 7 emissions. This fact has been used 
to propose a signal for the observation of v x neutrinos in water Cerenkov detectors [13]. Schematically the detection 
scheme works as follows (see Fig. 6). Supernova v x neutrinos, with average energies of 25 MeV, will predominantly 
excite 1~ and 2~ giant resonances in 16 via the le O(v x , v' x ) 16 0* neutral current reaction [10]. These resonances are 
above the particle thresholds and will mainly decay by proton and neutron emission. Although these decays will be 
dominantly to the ground states of 15 N and 15 0, respectively, some of them will go to excited states in these nuclei. 
In turn, if these excited states are below the particle thresholds in 15 N (E* < 10.2 MeV) or 15 (E* < 7.3 MeV), 
they will decay by 7 emission. As the first excited states in both of these mirror nuclei (E* = 5.27 MeV in 15 N and 
E* = 5.18 MeV in 15 0) are at energies larger than the SK detection threshold, all of the bound excited states in 15 N 
and 15 below will emit photons which can be observed in SK. 

Based on a calculation which combines the Continuum RPA with the statistical model [13], Supcrkamiokande is 
expected to observe about 700 7 events in the energy window E = 5 — 10 MeV, induced by v x neutrinos (with a 
FD distribution of T = 8 MeV), for a supernova going off at 10 kpc (« 3 • 10 4 light years or the distance to the 
galactic center), see Fig. 7. This is to be compared with a smooth background of about 270 positron events from the 
v e + p — > n + e+ reaction in the same energy window. The number of events produced by supernova u x neutrinos via 
the scheme proposed here is larger than the total number of events expected from ^-electron scattering (about 80 
events [44]). More importantly, the 7 signal can be unambiguously identified from the observed spectrum in the SK 
detector, in contrast to the more difficult identification from ^-electron scattering. 

The cascades of decays, following the inelastic excitation of 16 O by atmospheric or supernova v x neutrinos, can also 
result in the production of /3-unstable nuclei. If the Qp values of these nuclei are above the observational threshold 
energy in SK (~ 5 MeV), these decays might be detectable, and since they are usually delayed, might offer an 
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additional characteristic signature of the neutrino-induced reactions. As possible candidates 16 N, 15 C, 12 B and 12 N 
have been identified [45,46]. 

Fig. 8 shows the partial (P e ,e + ) and (y e ,e~) cross sections leading to the /3-unstable 16 N, 15 C and 12 B and 16 F, 
15 and 12 N ground states in the final channel; these cross sections reflect the sum over all cross sections with 
particle-bound states of these nuclei as these excited states will fast decay to the ground state by 7 emission. Our 
calculations have been performed up to neutrino energies E v = 500 MeV. At higher energies the total cross sections 
can be obtained from a relativistic Fermi gas model [47], including, however, additional channels like pion production. 

The results make it possible to draw some interesting conclusions. While the total cross sections increase with 
neutrino energies, most of this increase goes into new channels which open up with increasing neutrino energy and the 
partial cross sections leading to definite states have the tendency to saturate. Thus our partial cross sections to these 
states, obtained for, say, E v — 500 MeV, can be used to derive upper limits for the corresponding branching ratios 
expected for atmospheric neutrinos which, on average, have even larger energies. The total and selected partial cross 
sections for the neutral- and charged-current reactions on 16 O for different supernova neutrino spectra are summarized 
in Tables II-IV. 

A detailed discussion of the various cross sections is given in [46]. In general, these RPA cross sections are signifi- 
cantly smaller than estimated in [45]. In particular, they indicate that, for typical atmospheric neutrino energies, the 
partial cross sections leading to /3-unstable nuclei is, unfortunately, too small so that the observation of the decay of 
these nuclei does not constitute an additional viable signal for neutrino-induced reactions in Supcrkamiokande. 

For supernova v e neutrinos, however, the 16 0(P e , e + ) 16 N reaction can produce an observable additional signal in 
Superkamiokande for supernovae from within our galaxy. The reaction leads to excited states in 16 N at rather low 
energies, which then dominantly decay by neutron emission. However, as can be seen in Table IV, a sizable fraction of 
the (i> e , e + ) reactions also excite particle-bound states in 16 N, followed then by the (3 decay of the 16 N ground state. 
Assuming the standard antineutrino supernova spectrum (with T = 5 MeV) the partial cross section of 3.5 • 10~ 43 
cm 2 corresponds to about 40 supernova ^-induced events in SK leading to the 16 N ground state and which can 
be identified by the delayed (3 decay for a hypothetical supernova in the galactic center. Note, however, that this 
event rate corresponds to less than 1% of the total supernova neutrino rate in SK, with positrons being produced by 
the v e + p — > e+ + n reaction giving the dominating signal. Thus, it is unlikely that (3 decays from 16 N, generated 
by 16 0(z? e , e + ) 16 N in the supernova 1987A, were observed by the Kamiokande detector. The partial 16 0(P e , e + ) 16 N 
reaction cross section increases by more than a factor of 6, if complete v e <-> z^ T oscillations occur, constituting then 
a very sizable and clean signal for SK. 

Another interesting reaction, leading to an observable (3 decay in SK, is 16 0(P e , e + a) 12 B. For the standard spectrum 
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FIG. 7. Signal expected from supernova neutrinos in a water Cerenkov detector calculated for two different types of neutrino 
distributions (without chemical potential (above) and with chemical potential a — 3T and temperatures T = 6.26 MeV (for 
v x ) and T = 4 MeV (for v e )). The bulk of the signal stems from v e neutrinos reacting with protons, while the v x neutrinos 
induce the superimposed signal at energies E = 5 — 10 MeV (from [13]). 



we find a partial cross section of 1.1 • 1CT 44 cm 2 , which increases to 1.9 • 10~ 43 cm 2 for the case of oscillations. This 
translates into O(10) 12 B decays in the SK detector for a supernova in the galactic center at 10 kpc and if complete 
v e ^h,t oscillations occur. 



B. Argon 

The proposed ICARUS detector uses a liquid Argon technique providing energy and direction measurements for 
electrons, muons, pions and protons [48]. Observation of supernova neutrinos is also anticipated in the detector. 

The detection and analysis of supernova neutrinos by the ICARUS detector requires the knowledge of the neutrino- 
induced cross sections on 40 Ar for neutrinos and antineutrinos with energies up to about 100 MeV. At low neutrino 
energies the (v e ,e~) reaction is dominated by allowed GT transitions. However, the calculation of this cross section 
constitutes a quite challenging nuclear structure problem. In the Independent Particle Model (IPM), 40 Ar (Z = 18, 
iV = 22) corresponds to a 2-hole proton configuration in the sd shell and a 2-particle neutron configuration in the pf 
shell. One thus has to expect that cross-shell correlations will strongly influence the structure of the low-lying states 
in 40 Ar, including the GT_ response of the ground state. To describe such correlations requires shell model studies 
in the complete sd — pf shell which are currently not feasible by diagonalization methods. One therefore has to rely 
on truncated shell model calculations in which only a restricted number of particles are allowed to be excited from 
the sd shell to the pf shell. The best study to date has been performed by Ormand [49] who calculated the low-lying 
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FIG. 8. Total and partial cross sections to selected /3-unstable nuclear ground states for (v e ,e~) (upper), (i> e ,e + 
and (v, v') (lower part) reactions on le O as function of neutrino energy. As explained in the text, the symbols 15 0* 
in the lower panel denote the partial cross sections leading through particle-bound excited states in these nuclei. 
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GT response of the 40 Ar ground state within a 2hu shell model diagonalization in which he considered 2-particle 
excitations from the sd into the pf shell, carefully avoiding spurious center-of-mass excitations. His calculation is 
appropriate for solar neutrinos with energies E v < 14 MeV. Importantly, the shell model study showed that the 
GT transitions dominate over the allowed Fermi transition to the IAS in 40 K at 4.38 MeV. However, due to the 
(rather severe) truncation of the model space this shell model calculation violates the Ikeda sum rule and hence 
misses GT strength at higher excitation energies in 40 K. This strength will be important for neutrinos (like those 
from a supernova) which have larger energies than solar neutrinos and can excite the daughter nucleus at higher 
energies. An appropriate shell model calculation which describes this GT strength at higher energies is currently 
not available. However, RPA studies of the 4Q Ar(i/ e , e~) 40 K reaction have been performed considering allowed and 
forbidden multipoles up to J = 4. The respective cross sections are shown in Fig. 9. It is evident that GT transitions 
dominate the (y e ,e~) cross sections for neutrino energies E u < 50 MeV; at higher energies forbidden (in particular 
spin-dipole) transitions cannot be neglected. 

In the IPM, the GT+ strength for 40 Ar vanishes as all GT transitions, in which a proton is changed into a neutron, 
are Pauli-blocked. Although cross-shell correlations might introduce a non-vanishing GT+ strength, it should be 
small. Hence it is reasonable to calculate the (i> e ,e + ) cross section on 40 Ar within the RPA approach. The obtained 
results are also shown in Fig. 9. 



C. Iron 



The KARMEN collaboration used its sensitivity to the 56 Fe(z/ e , e~) 56 Co background events to determine the cross 
section for this reaction for the DAR neutrino spectrum and obtained a — (2.56 ± 1.08 ± 0.43) • 10~ 40 cm 2 [50]. We 
calculate a result in close agreement, a — 2.4 - 10 ~ 40 cm 2 [42]. 

We extended this investigation to the study of the charged- and neutral current reactions on 56 Fe. To allow also 
for the exploration of potential oscillation scenarios we also evaluated the cross sections and the knockout neutron 
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TABLE II. Partial cross sections for neutral-current neutrino-induced reactions on 1> 0. A Fermi-Dirac distribution with 
T = 8 MeV and zero chemical potential, which is typical for supernova and v T neutrinos and their antiparticles, has been 
assumed. The cross section, in units of lCP 42 cm 2 , represents the average for neutrino and antincutrino reactions, and the 
exponents are given in parentheses. The asterisks indicate that the cross sections have been summed over all particle-bound 
states. 
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TABLE III. Partial cross sections for charged-current neutrino- induced reactions on O. Fermi-Dirac distributions with 
T — 4 MeV and T = 8 MeV and zero chemical potential have been assumed. The first is typical for supernova v e neutrinos, 
while the second can occur for complete v e <-> oscillations. The cross sections are given in units of 10 _42 cm 2 , exponents are 
given in parentheses. 
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yields for various supernova neutrino spectra. Table V summarizes the total and partial cross sections for neutral 
current reactions on 56 Fe. For 56 Fe the neutron and proton thresholds open at 11.2 MeV and 10.18 MeV, respectively. 
But despite the slightly higher threshold energy, the additional Coulomb repulsion in the proton threshold makes 
the neutron channel the dominating decay mode. The total and partial cross sections for charged current (v e ,e~) 
and (^ e ,e + ) reactions on 56 Fe are listed in Table VI. For the standard supernova v e spectrum ((T, a) = (4,0)), the 
low-energy excitation spectrum is relatively strongly weighted by phase space. Hence, in that case the f e -induced 
reaction on 56 Fe leads dominantly to particle-bound states (~ 60%) and therefore decays by 7 emission. 

The work of Rcf. [42] has been extended to other iron isotopes, 52 ~ 60 Fe, by Toivancn et al. [51]. The respective total 
charged-current (u,e~) and neutral-current (v,v') cross sections for typical supernova neutrino spectra (i.e. T = 4 
MeV for v e neutrinos and T = 8 MeV for v x neutrinos) are shown in Figs. 10 and 11. Additionally these figures 
show the partial cross sections for the decay into the proton and neutron channels. Since for the isotopes 52 ~ 56 Co the 
proton threshold is lower than the neutron threshold, the preferred decay mode for the (v e , e~) reaction on 52 ~ 56 Fe is 
by proton emission. This is reversed for 58_60 Co, where decay into the neutron channel is preferred. Because parts 
of the Gamow- Teller and the Fermi strengths are located below the particle thresholds, these states decay by gamma 
emission which accounts basically for the difference between the total cross section and the sum of the partial proton 
and neutron decay cross sections. Decays into the a channel are unfavored. 

For the neutral-current reactions decay by proton emission is favored in the proton-rich nuclei 52_54 Fe. However, 
as the neutron threshold decreases with increasing mass number along the isotope chain, while the proton threshold 
energy increases, the probability for decay into the neutron channel increases at the expense of decay by proton 
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TABLE IV. Partial cross sections for charged-current antineutrino-induced reactions on O. Fermi-Dirac distributions with 
T = 5 MeV and T = 8 MeV and zero chemical potential have been assumed. The first is typical for supernova v e neutrinos, 
while the second can occur for complete v e <-> oscillations. The cross sections are given in units of 10 _42 cm 2 , exponents are 
given in parentheses. 
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FIG. 9. Total (u e ,e~) (solid) and (i/ e ,e + ) (dashed) cross sections on 40 Ar, calculated within the RPA approach. The 
short-dashed line shows the forbidden contributions to the (u e ,e~) cross sections. 



emission. This is clearly reflected in the trend of the partial cross sections as function of neutron excess. Due to 
pairing the neutron threshold is higher in even-even nuclei than in odd-A, explaining the odd-even staggering in the 
partial neutron decay cross sections. Again, the difference of the total cross sections compared to the sum of partial 
proton and neutron decay cross sections gives the (v, 1/7) cross sections, caused mainly by the Gamow- Teller strength 
below the particle thresholds. 

We note that neutrino-induced reactions on nuclei in the iron mass region might also play a role during the collapse 
and the shock-revival phase of a supernova [52,53]. Under these conditions, the cross sections have to be evaluated 
at the finite temperature of the stellar environment involving excited states of the parent nucleus. Relevant cross 
sections are presented in [54,55]. 



D. Lead 



The Fermi and Ikeda sum rules both scale with neutron excess (N — Z). As the charged-current response induced 
by supernova v e neutrinos (with average energies around 12 MeV) are dominated by Fermi and GT transitions, the 
charged-current cross sections on lead is expected to be significantly larger than on other materials like iron or carbon. 
This makes 208 Pb an attractive target for a supernova neutrino detector. 

To calculate the relevant neutrino-induced cross sections on 208 Pb we note that convergent shell-model calculations 
of the GT strength distribution are not computationally feasible. Thus, unlike in 40 Ar and 56 Fe, the X* = 1 + response 
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FIG. 10. Total charged current cross section for v e on 52 ~ 60 Fe (filled circles connected by the full line) and the partial neutron 
(empty circles, dashed line) and proton (crosses, dotted line) spallation cross sections. The difference between the total and 
the two partial cross sections gives the gamma emission cross section since the a channel is negligible. The u e spectrum with 
T = 4 MeV and a = was used. 



has been evaluated within the RPA approach which fulfills the Fermi and Ikeda sum rules. As the Sp+ strength (in 
this direction a proton is changed into a neutron) is strongly suppressed for 208 Pb, the Ikeda sum rule fixes the Sp- 
strength. In the calculation the A" = 1+ strength in 208 Pb was renormalized by the universal quenching factor which, 
due to a very slight A-dcpcndcnce is recommended to be (0.7) 2 in 208 Pb [56]. Thus, the Ikeda sum rule reads in this 
case Sp_ — Sf3 + w S@_ = 3 • (0.7) 2 • (N — Z). For the other multipole operators no experimental evidence exists for 
such a rescaling and we have used the RPA response. The RPA calculation is described in details in Ref. [57], which 
also discusses the relevant neutrino-induced cross sections for 208 Pb assuming a pion-decay-at-rest spectrum. 

The total and partial cross sections for charged current (v e , e~) and (& e , e + ) reactions on 208 Pb are listed in Table 
IVD. The (f e ,e) cross section on 208 Pb is about 20 time larger than for 56 Fe. This is caused by the (N — Z) and 
by the strong Z-dependence of the Fermi function. Furthermore, as the IAS energy and the GT_ strength is above 
the neutron threshold in 208 Bi at 6.9 MeV, most of the (f e ,e) cross section leads to particle- unbound states and 
hence decays by the neutron emission. Table VIII summarizes the total and partial cross sections for neutral current 
reactions on 208 Pb. 

There have been other calculations of the neutrino-induced cross sections on 208 Pb. The first study, performed in 
[37] , has been critisized and improved in [58] . These authors estimated the allowed transitions to the charged-current 
and neutral-current cross sections empirically using data from (p,n) scattering and from the Ml response to fix the 
Gamow- Teller contributions to the cross section. Ref. [58] completed their cross section estimates by calculating the 
first-forbidden contributions on the basis of the Goldhaber- Teller model. They calculated cross sections which are 
somewhat larger than the RPA results of [57]. More results have been reported in [59-61] which agree well with the 
RPA results. The later approach [61] is particularly interesting as it uses the experimental GT distribution, which 
has recently been determined in Osaka [62] , and the peaks of the spin-dipole response to constrain the (Hartree-Fock 
+ Tamm-Dancoff ) calculation. Furthermore, the spreading and quenching of the GT response has been considered by 
coupling to 2p-2h configurations. Suzuki and Sagawa obtain (3.2 -10~ 39 cm 2 ) for the [v e , e~) cross sections, assuming 
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FIG. 11. Total neutral current cross section for v x on 52_60 Fe averaged between v and v (filled circles connected by the 
full line) and the partial neutron (empty circles, dashed line) and proton (crosses, dotted line) spallation cross sections. The 
difference between the total and the two partial cross sections gives the gamma emission cross section since the a channel is 
negligible. The v x spectra with T = 8 MeV and a = were used. 



a DAR neutrino spectrum, in close agreement of Kolbe's RPA result (3.62 -10~ 39 cm 2 ) [57]. 

Proposed detectors like LAND and OMNIS will detect the neutrons produced by the neutrino-induced reactions on 
208 Pb. An obvious, but not very sensitive neutrino signal is the total neutron count rate. The two detectors might 
also be capable of detecting the neutron energy spectrum following the decay of states in the daughter nucleus after 
excitation by charged- and neutral-current neutrino reactions. The expected spectra have been predicted in [57]. The 
neutron spectrum for the charged current reaction on 208 Pb is dominated by the Fermi transition to the IAS and by 
the GT_ transitions. To understand the neutron spectrum one has to consider the neutron threshold energies for 
one-neutron decay (6.9 MeV) and for two-neutron decay (14.98 MeV) in 208 Bi. Hence the IAS and the collective GT 
resonance (with an excitation energy of about 16 MeV) will decay dominantly by 2n emission, while the low-lying 
GT_ resonance at E x = 8 MeV decays by the emission of one neutron. This has significant consequences for the 
neutron spectrum. In the 2-neutron decay the available energy is shared between the two emitted particles, leading to 
a rather broad and structureless neutron energy distribution. As can be seen in Fig. 12 this broad structure is overlaid 
with a peak at neutron energy around E n = 1 MeV caused by the one-neutron decay of the lower GT_ transition. 
One expects that, due to fragmentation which is not properly described in the RPA calculation, the width of this 
peak might be broader than the 0.5 MeV-binning that has been assumed in Fig. 12. The relative height of the peak 
compared with the broad structure stemming from the 2n-emission is more pronounced for the (T, a) = (4, 0) neutrino 
distribution than for a potential (T, a) = (8, 0) v e spectrum as it might arise after complete v e <-> oscillations. 

Ideally, OMNIS and LAND should have the ability to detect potential neutrino oscillations. However, as has been 
shown in [58] , the total neutron rate is by itself not suitable to detect neutrino oscillations, even if results from various 
detectors with different material (hence different ratios of charged-to-neutral current cross sections, as discussed above) 
are combined. Ref. [58] points out that in the case of 208 Pb an attractive signal might emerge. Due to the fact that 
the IAS and large portions of the GT_ strength resides in 208 Bi just above the 2-neutron emission threshold, Fuller 
et al. discuss that the 2-neutron emission rate is both flavor-specific and very sensitive to the temperature of the 
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TABLE V. Total cross sections for neutral current neutrino scattering on Fe for different neutrino energy spectra rep- 
resented as Fermi-Dirac distributions. The cross sections are in units of 10 _42 cm 2 and are averaged over neutrinos and 
antineutrinos. 



(T,a) 


(4,0) 


(6,0) 


(8,0) 


(10,0) 


(3,-3) 


(4,-3) 


(6.26,-3) 


bb Fo{u, u'-y) bb Fe 


2.5 ( 0) 


9.8 ( 0) 


1.7 ( 1) 


2.8 ( 1) 


1.2 ( 0) 


4.4 ( 0) 


1.6 ( 1) 


56 Fe(^, i/n) 55 Fe 


8.9 (-1) 


6.7 ( 0) 


2.2 ( 1) 


5.0 ( 1) 


2.8 (-1) 


1.7 ( 0) 


1.4 ( 1) 


56 Fe(^, i/p) 55 Mn 


1.2 (-1) 


1.0 ( 0) 


3.6 ( 0) 


9.3 ( 0) 


3.4 (-2) 


2.3 (-1) 


2.2 ( 0) 


56 Fe(^, v'a) 52 Cr 


2.4 (-2) 


1.9 (-1) 


6.6 (-1) 


1.7 ( 0) 


6.4 (-3) 


4.4 (-2) 


4.0 (-1) 


56 Fe(//, i/)X 


3.6 ( 0) 


1.8 ( 1) 


4.3 ( 1) 


8.9 ( 1) 


1.5 ( 0) 


6.3 ( 0) 


3.3 ( 1) 



TABLE VI. Total cross sections for charged current neutrino scattering on Fe for different neutrino energy spectra repre- 
sented as Fermi-Dirac distributions. The cross sections are in units of 10~ 42 cm 2 . 



(T,a) 


(4,0) 


(6,0) 


(8,0) 


(10,0) 


(3,-3) 


(4,-3) 


(6.26,-3) 


b0 Fe(^,e-7) OD Co 


9.8 ( 0) 


3.1 ( 1) 


6.1 ( 1) 


1.3 ( 2) 


7.7 ( 0) 


2.1 ( 1) 


7.5 ( 1) 


56 Fe(^ e ,e-n) 55 Co 


7.5 (-1) 


8.0 ( 0) 


3.2 ( 1) 


8.1 ( 1) 


2.5 (-1) 


1.7 ( 0) 


2.0 ( 1) 


56 Fe(^ e ,e"p) 55 Fe 


5.4 ( 0) 


3.2 ( 1) 


9.7 ( 1) 


1.7 ( 2) 


9.2 (-1) 


5.1 ( 0) 


4.7 ( 1) 


56 Fe(z/ e ,e-a) 52 Mn 


6.1 (-2) 


9.7 (-1) 


4.8 ( 0) 


1.5 ( 1) 


3.0 (-2) 


2.1 (-1) 


2.9 ( 0) 


5S Fc(!/ e ,e-)X 


1.6 ( 1) 


7.2 ( 1) 


1.9 ( 2) 


4.0 ( 2) 


8.9 ( 0) 


2.8 ( 1) 


1.4 ( 2) 


! " , Fe(77 e ,e + 7 ) bt, Mn 


2.3 ( 0) 


8.4 ( 0) 


1.8 ( 1) 


3.1 ( 1) 


1.4 ( 0) 


4.1 ( 0) 


1.6 ( 1) 


56 Fe(77 e ,e+n) 55 Mn 


4.2 (-1) 


4.0 ( 0) 


1.6 ( 1) 


4.0 ( 1) 


1.2 (-1) 


7.9 (-1) 


9.1 ( 0) 


56 Fe(77 e ,e+p) 55 Cr 


4.5 (-3) 


6.2 (-2) 


3.2 (-1) 


9.9 (-1) 


9.3 (-4) 


8.5 (-3) 


1.5 (-1) 


56 Fe(77 e ,e+a) 52 V 


1.1 (-3) 


1.7 (-2) 


9.5 (-2) 


3.1 (-1) 


2.0 (-4) 


2.1 (-3) 


4.3 (-2) 


56 Fe(77 e ,e+)X 


2.8 ( 0) 


1.2 ( 1) 


3.4 ( 1) 


7.2 ( 1) 


1.6 ( 0) 


4.9 ( 0) 


2.5 ( 1) 



v e distribution. To quantify this argument, in Ref. [57] the cross sections for the 208 Pb(^ e , e~2n) 206 Bi reaction was 
calculated in a model combining the RPA for the neutrino-induced response with the statistical model for the decay 
of the daughter states. The partial cross sections of 55.7 x 10 -42 cm 2 for v e neutrinos with (T, a) = (4, 0) Fermi-Dirac 
distributions should increase significantly if neutrino oscillations occur, as pointed out by [58]. For example, one finds 
for total u e <-> oscillations partial 2n cross sections of 1560 x 10~ 42 cm 2 (for neutrino distributions with parameters 
(T, a) = (8, 0)). We remark that these numbers will be probably reduced, if correlations beyond the RPA are taken 
into account, as part of the GT_ distribution might be shifted below the 2n-threshold. 

One has to note that this 2-neutron signal will compete with the 2-neutron decay stemming from the neutral-current 
reaction and hence will reduce the flavor-sensitivity of the signal. Ref. [57] predicts that the combined 2n-signal 
resulting from neutral-current reactions for the 4 v x neutrino types is larger than the one from the charged-currcnt 
reactions. However, if neutrino oscillations occur the neutral-current signal is unaffected while the charged-current 
signal is drastically enhanced. This supports the suggestions of Ref. [58] that the 2n-signal for 208 Pb detectors might 
be an interesting neutrino oscillation signal. 

V. NEUTRINO NUCLEOSYNTHESIS 

When the flux of neutrinos generated by the cooling of the neutron star in a type II supernova passes through 
the overlying shells of heavy elements, substantial nuclear transmutations are induced, despite the extremely small 
neutrino-nucleus cross sections. Specific nuclei (e.g. 10,11 B, 15 N, 19 F) might be, by a large fraction, made by this 
neutrino nucleosynthesis [63,64]. These are the product of reaction sequences induced by neutral current {y,v'} 
reactions on very abundant nuclei such as 12 C, 16 or 20 Ne. If the inelastic excitation of these nuclei proceeds to 
particle-unbound levels, they will decay by emission of protons or neutrons, in this way contributing to nucleosynthesis. 
As the nucleon thresholds are relatively high, effectively only , v T neutrinos and their antiparticles with their higher 
average energies contribute to the neutrino nucleosynthesis of these elements. It has been noted that the neutrino 
induced nucleosynthesis, the so-called ^-process, might also be responsible for the production of 138 La and to a fraction 
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FIG. 12. Neutron energy spectrum produced by the charged-current (i/ e ,e~) reaction on 208 Pb. The calculation has been 
performed for different supernova neutrino spectra characterized by the parameters (T,a). Note that the cross sections for 
(T, a) = (4, 0) and (3,3) have been scaled by a factor of 5. 



of the 180 Ta abundance [63]. The 138 La nuclide is of special interest as it appears to be produced by the charged 
current (u e ,e~) reaction on the s-process element 138 Ba, as has already been proposed in [63,65] and recently been 
confirmed in detailed studies by Heger et al. [67]. This finding is quite welcome as it makes the v process sensitive 
to the flux and distribution of supernova v e and ^„,i/ T , complementing the constraints for supernova v e neutrinos 
from their observation in the water Cerenkov detectors for SN1987A. Neutrino nucleosynthesis is thus potentially an 
important test for the predictions of supernova models. This test is particularly stringent, if neutrino oscillations 
involving v e neutrinos occur in the supernova environment. 

Neutrino nucleosynthesis has been proposed by Woosley et al. [63] who also performed detailed production studies 
in a 20 M Q star. Later Timmes et al. have extended this calculation to studies of the v process in a full galactical 
model [64] . Very recently, Heger et al. [67] have improved these earlier studies by considering mass loss in the evolution 
of the progenitor stars, and by using a complete and updated reaction network which includes all the heavy elements 
through bismuth. This improvement is essential for a consistent treatment of the s-process in the progenitor star. 
Last, but not least the studies of [67] also used improved neutrino- induced reaction cross sections for the key nuclei of 
the v process. For the p and sd shell nuclei the allowed neutrino response has been evaluated on the basis of the shell 
model, while the forbidden responses were calculated within the framework of the RPA. For the heavier nuclei, in 
particular for the progenitors of 138 La and 180 Ta, all cross sections were derived within the RPA. The various partial 
decay cross sections were then evaluated using a statistical model cascade. 

As stressed above, Heger et al. [67] followed the neutrino nucleosynthesis of the light (e.g. n B, 19 F) and heavy 
(e.g, 138 La, 180 Ta) candidate nuclei in a self-consistent way in complete stellar evolution models that included the 
evolution of all isotopes up to bismuth from the time the star ignited central hydrogen burning through the supernova 
explosion. A consistent stellar modelling, including a sufficient nuclear network, is important for two reasons. At 
first, it is essential to reliably describe the production of the progenitor nuclides (e.g. 12 C, 20 Ne, 138 Ba...) in the 
presupernova evolution. Secondly, the nuclides produced by neutrinos from these progenitors can suffer severe photo- 
dissociations due to heating by the supernova shock wave which passes the region of neutrino nucleosynthesis after 



20 



TABLE VII. Total cross sections for charged current neutrino scattering on nuclei for different neutrino energy spectra 
represented as Fermi-Dirac distributions. The cross sections are given in units of 10~ 42 cm 2 . 



(T,a) 


(4,0) 


(6,0) 


(8,0) 


(10,0) 


(3,3) 


/ A 0\ 

(4,3) 


(6.26,3) 


m Pb(^ e , e - 7 f J8 Bi 


4.7 ( 1) 


1.3 ( 2) 


2.5 ( 2) 


4.0 ( 2) 


3.5 ( 1) 


7.6 ( 1) 


2.2 ( 2) 


208 Pb(^,e-n) 207 Bi 


2.3 ( 2) 


9.9 ( 2) 


2.3 ( 3) 


4.0 ( 3) 


1.2 ( 2) 


4.2 ( 2) 


1.9 ( 3) 


2U8 Pb(^ e ,e p) 2U 'Pb 


1.8 (-2) 


1.1 (-1) 


3.3 (-1) 


6.9 (-1) 


7.2 (-3) 


3.3 (-2) 


2.3 (-1) 


208 Pb( I , e ,e-a) 204 Tl 


2.1 (-2) 


2.6 (-1) 


1.1 ( 0) 


3.0 ( 0) 


4.7 (-3) 


4.1 (-2) 


6.0 (-1) 


208 PbK,e")X 


2.8 ( 2) 


1.1 ( 3) 


2.5 ( 3) 


4.5 ( 3) 


1.6 ( 2) 


4.9 ( 2) 


2.1 ( 3) 


2U8 Pb(77 e , e + 7 ) 2,J8 Tl 


5.8 (-1) 


3.0 ( 0) 


7.9 ( 0) 


1.5 ( 1) 


2.7 (-1) 


1.1 ( 0) 


6.1 ( 0) 


208 Pb(77 e ,e + n) 207 Tl 


4.9 (-1) 


3.8 ( 0) 


1.5 ( 1) 


3.9 ( 1) 


2.0 (-1) 


8.9 (-1) 


8.5 ( 0) 


208 Pb(F e ,e+p) 207 Hg 


1.7 (-7) 


1.4 (-5) 


2.2 (-4) 


1.5 (-3) 


8.4 (-9) 


3.2 (-7) 


4.2 (-5) 


208 Pb(77 e ,e+ a ) 204 Au 


4.3 (-7) 


4.0 (-5) 


6.5 (-4) 


4.4 (-3) 


2.1 (-8) 


8.1 (-7) 


1.2 (-4) 


208 Pb(F e ,e+)X 


1.1 ( 0) 


6.8 ( 0) 


2.3 ( 1) 


5.4 ( 1) 


4.7 (-1) 


1.9 ( 0) 


1.5 ( 1) 



TABLE VIII. Total cross sections for neutral current neutrino scattering on Pb for different neutrino energy spectra 
represented as Fermi-Dirac distributions. The cross sections are in units of 10 _42 cm 2 and are averaged over neutrinos and 
antineutrinos. 



(T,q) 


(4,0) 


(6,0) 


(8,0) 


(10,0) 


(3,-3) 


(4,-3) 


(6.26,-3) 


2U8 Pb(^' 7 ) 2U8 Pb 


1.4 ( 0) 


7.4 ( 0) 


2.1 ( 1) 


4.5 ( 1) 


7.0 (-1) 


2.5 ( 0) 


1.5 ( 1) 


208 Pb(^'n) 207 Pb 


1.2 ( 1) 


4.8 ( 1) 


1.2 ( 2) 


2.3 ( 2) 


6.9 ( 0) 


2.0 ( 1) 


9.4 ( 1) 


208 Pb(^'p) 207 Tl 


1.6 (-5) 


3.5 (-4) 


2.4 (-3) 


8.7 (-3) 


2.9 (-6) 


3.1 (-5) 


9.0 (-4) 


208 Pb(^'a) 204 Hg 


7.8 (-5) 


3.0 (-3) 


2.6 (-2) 


1.1 (-1) 


8.1 (-6) 


1.5 (-4) 


7.9 (-3) 


208 Pb(^X)X 


1.3 ( 1) 


5.6 ( 1) 


1.4 ( 2) 


2.7 ( 2) 


7.6 ( 0) 


2.3 ( 1) 


1.1 ( 2) 



the neutrinos. The neutrino production factors for 4 key nuclides are summarized in Fig. 14 for 15 Mq and 25 Mq 
stars. The results are shown relative to the 16 production in these stars recognizing that the solar 16 abundance 
is an indicator for core-collapse supernova element production. 

The calculation [67] confirms earlier studies [63,64] that the neutrino nucleosynthesis makes a large fraction of the 
solar 11 B and 19 F abundance. For both nuclides, neutral-current reactions, on 12 C and 20 Ne, induced by v x neutrinos 
are the main production source. The new calculation predicts a somewhat smaller production of 19 F than the earlier 
studies. A measurement of the GT response on 20 Ne might be quite desirable, including the cascade of decays of 
excited states, to better constrain the (v, v') cross section on 20 Ne. As predicted in [63,65], 138 La is mainly produced 
by charged-current reactions on 138 Ba, while the 7-process contribution is small and the neutral-current contribution 
from 139 La, which had been speculated to be the possible 138 La production process via (v, v'n), is insignificant. One 
observes that enough 138 La is being made to explain the solar abundance, where the s-process production of the 
parent 138 Ba in the s-process prior to the supernova plays an essential role. The calculations also show a significant 
production of the rarest nuclide 180 Ta, where both, charged- and neutral-current processes contribute. The 180 Ta 
production factors in Fig. 14 adds the contributions of the ground state and the 9~ isomeric state. As the ground 
state is rather short-lived (~ 4.5 h), the isomeric state is the sole contributor to the solar 180 Ta abundance. Estimates 
[66] indicate that at thermal freeze-out about 30 — 50% of 180 Ta is in the isomeric state bringing the production factors 
of Fig. 14 in closer agreement with the solar abundances. 

Importantly the 138 La and 180 Ta production is sensitive to the neutrino distributions. For example, if the tempera- 
ture of the v e spectrum increases from 4 MeV to 6 MeV (which might correspond to a neutrino oscillation scenario) the 
138 La production factor increases by a factor 2 and, in the study of [67], this nuclide would already be overproduced. 
This makes the neutrino nucleosynthesis of 138 La a potentially interesting test for neutrino oscillations, in particular 
for the yet unmeasured neutrino mixing angle 9\% [67]. 

Finally, we note that v nucleosynthesis of 138 La and 180 Ta competes with the p-process production of these elements, 
making a reliable determination of the p-process abundances also important. 
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FIG. 13. Partial cross sections for the ^-induced charged-current reaction on 138 Ba. The calculations were performed for 
Fermi-Dirac neutrino spectra with a = and different temperature values (from [67]). 



VI. CONCLUSIONS 

In this review we have described the methods used in the evaluation of the neutrino- nucleus reaction cross sections. 
We have shown, in detail, the results of calculations of these cross sections for a variety of targets of practical 
importance. 

Many of the described reactions are accessible in experiments with a neutrino source from the pion and muon 
decay at rest, available at the future very intense neutron spallation sources. Detailed comparison of the results of 
such experiments would establish important benchmarks for comparison of theory and experiment. This, in turn, 
would lead to refinements or possible modifications of the theoretical treatment of processes involving the charged 
and neutral current interaction of neutrinos with complex nuclei. Having a reliable tool for such calculation is of great 
importance in a variety of applications, like the study of neutrino oscillations, detection of supernova neutrino signal, 
description of the neutrino transport in supernovae, or description of the r-process nucleosynthesis. 
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(from [67]). The open (filled) symbols represent stellar evolution studies in which neutrino reactions on nuclei were excluded 
(included). 
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